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Abstract. Motivation is given for trying a theory of gravity with a preferred ref- 
erence frame ("ether" for short). One such theory is summarized, that is a scalar 
bimetric theory. Dynamics is governed by an extension of Newton's second law. In 
the static case, geodesic motion is recovered together with Newton's attraction field. 
In the static spherical case, Schwarzschild's metric is got. An asymptotic scheme of 
post-Minkowskian (PM) approximation is built by associating a conceptual family of 
systems with the given weakly-gravitating system. It is more general than the post- 
Newtonian scheme in that the velocity may be comparable with c. This allows to justify 
why the 0PM approximation of the energy rate may be equated to the rate of the New- 
tonian energy, as is usually done. At the 0PM approximation of this theory, an isolated 
system loses energy by quadrupole radiation, without any monopole or dipole term. 
It seems plausible that the observations on binary pulsars (the pulse data) could be 
nicely fitted with a timing model based on this theory. 



1 Introduction 

The dominant opinion is that we already have both an excellent theory of grav- 
itation - namely, Einstein's general relativity (GR), and an excellent theoretical 
framework for particle physics - namely, the standard model, which is based on 
quantum field theory (QFT). The important problem of making quantum the- 
ory and GR compatible together should then be solved essentially by going to 
a superstructure including the two former theories as particular and/or limiting 
cases, the candidates for that superstructure being searched in the wide range 
of string theory. This might prove to be the case in the future, but it must be 
admitted that string theories are very complex, being based on differential geom- 
etry on high-dimensional manifolds. Even GR and QFT separately are already 
very complex, to the point that the link between either theory and its exper- 
imentally confirmed predictions is not fully understood, being based on rather 
peculiar algorithms Therefore, apart from the main line of research, one can 
find some risky attempts, in which it is tried to take a quite different route. 



One such route is suggested by the existence of an alternative version/ in- 
terpretation of special relativity, which version was initiated by Lorentz ^2^ and 
Poincare |3l4j . According to this version, the "relativistic" effects, which essen- 
tially follow from the Lorentz transformation, are all due to the "true" Lorentz 
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contraction of physical objects as they are moving through the "ether" or fun- 
damental inertial frame [515] , In fact it can be shown that both Lorentz's space- 
contraction and Larmor's time-dilation follow from the negative result of the 
Michelson-Morley experiment, under the assumption that light propagate in all 
directions with the same velocity c with respect to a rigid reference frame, in 
which Euclidean geometry applies |6I7) . In this "Lorentz-Poincare version of 
special relativity", the time which is defined in an inertial frame by using the 
Poincare-Einstein clock synchronization is generally not the "true time" : only 
in the fundamental inertial frame does the Poincare- Einstein-synchronized time 
coincide with the true time |8l9l6j . The Lorentz-Poincare version of special rel- 
ativity has the same equations as the usual version initiated by Einstein in 
his well-known June 1905 paper, and is observationally equivalent to the usual 
version However, just because of this, the "ether", which is still there in 
the former version, is indetectable. In particular, the true time of the Lorentz- 
Poincare version is not experimentally accessible. 

But, precisely, an essential assumption of the investigated alternative theory 
of gravitation |10l7llllf^ is that the presence of a gravitational field breaks the 
Lorentz symmetry, thus making the ether detectable Therefore, in this the- 
ory, the inertial time in the ether frame is a physically identifiable preferred time 
T (called the "absolute time"). It is accepted that the existence of a preferred 
space-time foliation could give a solution to extend quantum theory to the situa- 
tion with gravitation [1'6\ . The very beginning of this program is the formulation 
of the free Klein-Gordon wave equation in a gravitational field, which is ambigu- 
ous in GR, but leads to a unique equation (which is a preferred-frame one) in 
the framework of the "ether-theory" |14j . Another motivation for trying a very 
alternative theory of gravitation is to solve some problems which are common to 
GR and to most extensions of it - namely, the existence of singularities (there 
is indeed no singularity in the investigated theory, neither during gravitational 
collapse "in free fall" with spherical symmetry |15| nor in homogeneous cosmo- 
logical models ^E]), and the interpretation of the gauge condition (there is no 
gauge condition in this scalar theory). 

This theory gives the same predictions as GR for light rays jJT] and it un- 
ambiguously predicts that the cosmic expansion must be accelerated ^B]. As 
to celestial mechanics: since this theory predicts Schwarzschild's motion for test 
particles in the static case with spherical symmetry 'T'll^, one a priori expects 
that it should improve celestial mechanics as compared with Newton's theory, 
in the same way as GR. However, due to the existence of preferred-frame ef- 
fects, one has to check this carefully. To this end, an "asymptotic" scheme of 
post-Newtonian (PN) approximation has been built for a weakly self-gravitating 
system of extended bodies, that consists in associating a one-parameter family 
of systems with that given system: each system is defined by an initial-value 
problem [18; . This is a correct way to obtain asymptotic expansions of the solu- 
tion fields, with respect to the small weak-field parameter. For extended bodies 
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in GR, this has been done only with a very particular initial condition as to the 
space metric jj^li and the obtained local PN equations have not been integrated 
in the volume of the bodies to get equations of motion of the mass centers. The 
latter has been done for the investigated scalar theory |2()I21| . The obtained 
equations of motion of the mass centers depend on the spin rates of the bodies 
and involve structure-dependent parameters - two natural features for a theory 
that integrates the mass-energy equivalence. In order to test the theory, these 
equations have been implemented 1221281 in a code for numerical integration 
and parameter adjustment, that has been built by the author and tested t24u2Eij 
with standard equations of motion. (For the moment, the spin rates of the main 
bodies of the solar system have been neglected.) It has thus been found that 
the difference with a standard ephemeris does not exceed 3" over the 20th 
century [53]. One may expect to significantly improve the observational agree- 
ment by improving the numerical treatment and also by adjusting the equations, 
not on an ephemeris, but directly on observational data. Indeed some correction 
factors of observational data are taken as free parameters in the adjustment 
of an ephemeris, hence the observations are not completely independent of the 
gravitational model. 

Thus the current situation as to the experimental test of this theory is en- 
couraging. However, it is usually considered that a scalar theory should have 
problems with gravitational radiation (Note that it is usually considered 
that a scalar theory should also have problems with gravitational effects on light 
rays |2H1, and that it is not the case for the present preferred-frame theory.) 
Therefore, the aim of the present investigation was to examine gravitational ra- 
diation. Before discussing the main line of that investigation, we shall present 
the theory. 

2 Basic Principles of this Scalar Theory 
2.1 Space-time metric and field equation 

The theory considered is a bimetric theory: space-time is endowed with both a 
flat metric and a curved metric 7. The relation between the two metrics is 
flxed by just a scalar field / j7j. This can be true only because this is a preferred- 
frame theory and the equations are written in the preferred reference frame E 
("ether"). ^ Thus, space-time is the product R x M, where M is the "space" 
manifold, i.e. the set of the positions x in the preferred reference frame E, and 
most equations are merely space-covariant. The preferred time coordinate T is 

^ The theory could be rewritten in a generally-covariant form by incorporating the 
velocity field of the ether (relative to the reference frame associated with the given 
coordinate system) into the fields of the theory. If one restricts the choice of the 
reference frame to the "inertial" class made of the frames that have a uniform rigid 
motion with respect to E (the uniformity and rigidity being defined in terms of the 
fiat metric 7"), then one just has to add a constant vector V to the unknowns. 
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the inertial time in the frame E (i.e. T — /c with c the velocity of hght, and 
where the coordinates {x^){iJi = 0, ...,3) are Gahlean for the metric 7*^, that 
is, I'^u = Vfj.i' with (yy^jy) = diag(l, — 1, — 1, — 1)). In particular, the scalar field 
/ = (700 )e (thus in coordinates (j/'^) bound to the frame E and with y° = cT) 
determines the slowing down of real clocks in a gravitational field, as compared 
with an "absolute clock" that would measure the absolute time T. This is the 
case for a clock that is fixed in the ether frame, and far enough from any massive 
body so that it is not affected by any gravitational field [J] . 

The equation for the field / is a nonlinear wave equation: 



with A = divgogradgo the usual Laplace operator defined with the Euclidean 
metric g", which is the spatial part of the flat metric 7° in the frame E; G is 
Newton's gravitation constant and a = (T*'°)e is the mass-energy density in the 
ether frame, T being the "mass tensor" , i.e. the energy-momentum tensor in 
mass units. In the most general case, applicable to a heterogeneous universe on 
a cosmological time-scale, the gravitational field is not directly the field / but 
instead the field of the "ether pressure" Pe, and a different field equation applies 
[7]. But this reduces to (Q) if the time variation of the "reference pressure" 
p^{T) = SupxgM?'e(x, T) is neglected ^TD, which should be the case except on a 
cosmological time-scale. The analysis of homogeneous cosmological models (16) 
confirms that the variation of p^(T) takes place over long time scales, of the 
order of 10* years. 

2.2 Dynamical equations 

Motion of a free test particle is defined by an extension of the special-relativistic 
form of Newton's second law: 



where E is the energy of the test particle: for a mass particle, E = m{v)c^ with 
m{v) — 7^m(0) the velocity-dependent inertial mass (71, is the Lorentz factor); 
for a photon, E = hv with v the frequency measured with real clocks affected 
by the gravitational slowing down and thus measuring the "local time" ix (see 
Eq. Q below). Moreover, 





(1) 



c 



(3) 



is the gravity acceleration, 



dx 1 dx 



(4) 



V = 



dix VJdT 



Gravitational Energy Loss in a Scalar Theory of Gravity 



5 



is the velocity, and D /Dt-^ is the correct time-derivative of a vector in the space 
manifold M endowed with the time-dependent Riemannian metric g (g is the 
spatial part of 7 in the frame E) Finally, 

P-fv (5) 

is the momentum of the test particle. 

In the static case {f^ — 0), the extension ||2Jl of Newton's second law (thus 
three scalar equations) implies that free test particles follow the geodesic lines of 
metric 7 (thus four scalar equations: the time component of the geodesic equa- 
tion is equivalent to the energy equation valid for a test particle in a static field) 

In the static case with spherical symmetry, Schwarzschild's exterior metric 
is obtained [Jj. 

The dynamical equation for a continuous medium (fluid, electromagnetic 
field, ...) is deduced from Newton's second law as defined above for a test particle. 
Indeed, for a dust, Eq. (0 may be applied pointwise and implies [T^ : 

T;.,. = b,, (6) 

&o(T)EEi5,,.,oT^"^ 5.(T)^-ig,fc,oT'"= (7) 

(recall that g = (gij) is the (curved) space metric in the frame E. Semicolon 
means, as usual, covariant differentiation defined with the connection associated 
to the (curved) space-time metric 7. Indices are raised and lowered with 7). 
The universality of gravity means that the same equation must hold true for any 
continuous medium. 



2.3 General energy conservation in this theory 

Equation © may be rewritten in terms of the flat metric 7^^, giving as the time 
component: 

ro^.. = ^/,oT"° (8) 

(valid in Galilean coordinates for the flat metric) ^21- Using the equation ^ 
for the scalar field /, the r.h.s. of Eq. © may be transformed to a 4-divergence 
(with respect to the flat metric), thus giving a true conservation equation for 
the energy density 

= £m + Eg- (9) 

The energy rate in a fixed domain SI, whose boundary dil is not crossed 
by any matter, is got by integrating the conservation equation, it is 

E = -^, f=-4^/ drfg.ndS, E= [ edV. (10) 



SttG 
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3 Gravitational Radiation 

In GR and other relativistic theories of gravitation, gravitational radiation is 
mostly studied with the "weak-field approximation" , which is essentially a stan- 
dard linearization of the field equations j29l8()i:ilt^ . However, we are here in 
a domain where the accuracy of the observations (of the timing of the pulses 
emitted by binary pulsars) is very high, and so also is the precision of the com- 
parison theory/ observation |32]- Therefore, it is desirable that the link between 
the equations of the theory and its numerical predictions should be obtained 
with the help of a well-defined approximation scheme. Thus one would need 
to introduce an asymptotic framework, associating with the physical system of 
interest S a family (Sa) of gravitating systems (A is the small field-strength pa- 
rameter). A such "asymptotic scheme" was previously developed for celestial 
mechanics (see the fourth paragraph of the introduction here). As A ^ 0, the 
behaviour becomes Newtonian. Thus it is a post- Newtonian (PN) scheme. Like 
the standard PN scheme |l29l38l81IH4j . that "asymptotic" PN scheme leads to 
Poisson equations with instantaneous propagation, therefore it does not fit in the 
present context, since the aim is to describe how the energy balance is affected 
by the propagation of gravitational waves at the finite velocity c. 

What one would need is an asymptotic post-Minkowskian (PM) scheme, i.e., 
one leading to a wave equation for the gravitational field (here the scalar /). 
In GR, it has been proved that one-parameter families of solutions to the vac- 
uum Einstein equations generically exist, and that the successive coefficients of 
their Taylor expansions with respect to the parameter do satisfy the successive 
equations of the PM approximation [35]. However, in order to justify standard 
calculations based on equating the Newtonian energy rate to the energy rate as 
given by the "quadrupole formula of GR" (see e.g. Refs. 30 31 34j), it is neces- 
sary that there is indeed matter which has the relevant quadrupole tensor. One 
possible approach that would provide a reasonable justification to such calcu- 
lations would be to show that the same given self-gravitating system may be 
consistently envisaged both in well-defined PN and PM approximation schemes; 
that the PM scheme (as applied to the given system and thus also inside the 
bodies) gives some calculable non-zero gravitational energy rate EpM', and that 
the PM scheme is more general than the PN scheme, so that it does make sense 
to state the equation 

E ~ EpM — -EpN- (11) 

This turns out to be feasible. We shall insist on the essential points of that 
approach (the details of which shall be published elsewhere) and shall skip the 
calculations. 
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3.1 Principle of the asymptotic PM approximation 

As for the "asymptotic" PN scheme, which has been presented in detail in Rcf. 
| 18) . a family (Sa) of perfect-fluid systems is defined by a family of initial condi- 
tions : at the initial time, the fields of pressure and proper rest-mass density in 
system Sa have respectively the form 

p(^)(x)-AV'Hx), p(^)(x) = Ap(i)(x); (12) 

the initial condition for the gravitational field /, or rather for V = (c^/2)(l — /) 
(which plays the role of the Newtonian potential) , is given by 

y(^)(x) = AV^(i)(x), dTV^^\^) = XdTV'^^\x); (13) 

and the initial condition for the velocity is 

u(^)(x) = u(^)(x). (14) 



The system of interest, S (e.g. a binary pulsar), is assumed to correspond to 
a small value Aq <C 1 of the field-strength parameter 

A = Sup^^mII - /^^Hx,r = 0)]/2 = Sup,eM[^^^Hx,T = 0)]/c\ (15) 

Hence, equations C2Ji-2 and H13|l i mean physically that, in the system S (and 
more precisely near the center of the most massive body of S), the scalar 1 — / = 
1 ~ f^'^°\ that determines the magnitude of the difference between the flat and 
curved metrics, is of the order of magnitude 2Ao, and that the fields p and p have 
the order Ag and Aq, respectively (in units such that the fields p'^^-' and p'^^^ are 
of the order of unity near the center of the most massive body). These equations 
()12|li-9 and 1)13(1 1 are exactly the same as in the PN scheme JHl- ^ The essential 
difference with the (asymptotic) PN scheme of the scalar theory JH| is that, in 
the PM scheme, it is not necessary that |u| <C c (in the physical sense) for the 
given system S, because the velocity in the generic system Sa of the family is 
order zero with respect to the parameter A (Eq. 1(14(1 ') - whereas the velocity 
is order \/A in the PN scheme Correspondingly, the time derivative of the 
gravitational potential V is order 1 in A in the PM scheme (Eq. ((13(l 9^. whereas 
it is order 3/2 in the PN scheme ^Sl- Hence, the (asymptotic) PM scheme is 
just more general than the corresponding PN scheme, in the sense that the given 
system S will automatically be suitable for the PM scheme if it is for the PN 
one, and the converse is not true. (But the PM family (Sa) is not "more general" 
than the corresponding PN family: it is just a different family of solutions.) Thus, 
EpM can be used in the PN scheme. 



^ The small parameter considered in the present paper is A = , where e is that used 
in Ref. 
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3.2 PM approximation of the energy rate 

We assume that the initial-value problem Pa, defined by Eqs. H12|l -114 () . has 
generically a unique solution in some large-enough domain, and that the solu- 
tion fields admit asymptotic expansions in powers of A, which have necessarily 
the same dominant order in A as the initial values. Precisely, we need just the 
principal part in these expansions. After changing the mass and the time units, 
for system Sa, to the new units [M]a = A[M] and [T]a = [T]/VA (where [M] 
and [T] are the starting units), the small parameter A becomes proportional to 
1/c^, and all fields: p, p, V, dV/dx^ , hut the velocity, are order zero, whereas the 
velocity field u is order c. We may hence write: 

V = U + 0{\/^), a = cro + 0(l/c2), u = uoc + 0(l/c) (16) 

(recall that a = (T°°)e is the source of the gravitational field, Eq. (^). These 
expansions have primarily to be valid at fixed values of the relevant time and 
space variables. Due to the fact that the velocity is order zero in the small 
parameter (in fixed units) for the PM family, the relevant time variable (the 
"dynamical time" in the terminology of Ref. 19 ) is the true time T (as expressed 
in fixed units) - whereas the relevant time variable is T\/\ for the PN family 
of systems |18) . Therefore, in the new (varying) units, in which the expansions 
have the convenient form (|16|l . and in which A (x 1/c^, the relevant time variable 
is = cT. As a consequence, the expansion of the field equation leads, as 
we wished, to the wave equation for the 0PM potential U: 

UU = - AU = AnGaa (x" = cT) (17) 

- instead of the Poisson equation which applies to the OPN (Newtonian) poten- 
tial. As A ^ 0, the gravitational energy fiux H1U|I has a principal part 'Pqpm- 

<? = <?opm(1 + 0(A)), <PoPM = ^f drUUku'^dS, (18) 
where U is the relevant (retarded) solution of Eq. (fTTjl . 



3.3 There are only quadrupole terms in the Hmit radiative flux 

The next task is to calculate the limit, as i? ^ cx), of the 0PM flux «?opm on the 
sphere dSl : |X| = R. To this end, one first expands the retarded potential U as 
R ^ oo, then one inserts this in H18() . The multipole expansion of the retarded 
potential U contains monopole, dipole, and quadrupole terms, i.e. it contains 

M=JcrodV, d,(X,r) EE /2;Vo(x,cT- |X|)dt/(x), (19) 
Jy (X,T) = / a;*x%o(x, cT - |X|)dF(x). (20) 



Indeed, we find: 



G R cR^ 2c2i?3 ••• \R2 
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where overdot means partial derivative with respect to T and the points of 
suspension indicate omitted terms of order coming from time derivatives 
of order > 3 in the Taylor series expansion of the retarded density. However, 
M — Const, disappears from the flux H18|l . due to the presence of the partial 
derivative drU . As to di, it is in the limit flux just by di. But = from the 
0PM conservation of the total momentum (which is obtained when inserting the 
expansions H16() into the dynamical equation ©). It follows that the limit as 
R ^ oo oi the PM energy flux contains just quadrupole terms: 

(^-OPHjiim = (2 tr P + tr2 j) . (22) 

3.4 Use of the PM approximation of the energy flux 

As already mentioned, the fact that the PM scheme is more general than the 
PN scheme implies that it makes sense to use the energy rate -EpM in the PN 
scheme. Let us try to expand on this point. The limit, as i? — > 00, of the principal 
part of the PM energy flux, is given by Ea. H22l) . whereas, at the zero-order PN 
approximation, the corresponding limit is zero. Indeed, the principal part of the 
PN energy flux, ^opNj is just the Newtonian energy flux. But the Newtonian 
energy of an isolated system is conserved, hence we have ^ 

hm <?>oPN = 0. (23) 

However, although the zero-order PN (i.e., Newtonian) energy -Bqpn is thus 
constant at this same zero-order PN approximation, it has no reason to remain 
so at the following, truly post-Newtonian, approximations. Indeed, as the order 

• (n) 

n of the PN approximation increases, the energy rate E^^ calculated at the 
j^th pj^ approximation should become a better and better approximation of the 
exact energy rate more precisely one should have 

i; = 4'^+0(A"+i). (24) 

' in) 

Moreover, the dominant part of Ep^ should precisely be the rate of i?oPN: if the 
PN expansions are uniform with respect to time, then, by summing the rates of 
the successive terms in the asymptotic expansion of the energy: 

E = £;oPN + A(J£;ipn) + ... + 0(A"+i), (25) 

one should obtain the asymptotic expansion of the energy rate: 

E - ^(£;opn) + A^(<5i?iPN) + ... + 0(A"+i). (26) 

^ This is easy to check by integrating the local balance equations of Newtonian gravity 
|11|: the point is that, in Newtonian gravity, the gravitational energy flux has just 
the same expression II8II 9 as in the 0PM approximation, but with U being now the 
Newtonian potential, which is such that U^k ~ 0{l/R^) and drU — ^ at large R. 
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More exactly, when we write Eq. H26I) . we are using the n*"^ PN approximation 
and we must hence rewrite this equation as 

i^ = 4"^+0(A"+i)= fA") (^Eopn) + x(^) (<5£;ipn) + ... + 0(A"+1). 

V"-' / nPN V"-' / nPN 

(27) 

Of course, just the same may be done with the PM scheme. What happens is 
that the first iterations of the PN scheme, in particular the very first one (the 
Newtonian approximation), give the simple prediction that the energy of the 
gravitating system is conserved, i.e. 

E^""^ =0 for < n < p. (28) 

Therefore, the exact energy rate is certainly very small, and one would like to 
calculate it by using the p PN approximation, that gives the first non-zero 
result (in GR, the standard PN scheme gives p = 5/2 36 , the fractional number 
meaning merely that a Taylor expansion in A does not apply any more) . But this 
is very complicated, hence one is led to admit that the first non-zero result given 
by the PM approximation, and which turns out to correspond to the zero-order 
PM approximation (Eq. (|22|l '). is a good approximation to the exact result. The 
reason to believe so is that the PM approximation, in contrast to the PN one, 
does take into account naturally the phenomenon of gravitational radiation. We 
may summarize by admitting that, for some number p (may be p = 5/2): 

E^^E^J'^c^f^) iEoP^)^f±) (Eopm). (29) 

V^-* / pPN V'^-' / 0PM 



3.5 Peters-Mathevirs coefficients and decrease in the orbital period 

Now we consider a binary system and we assume i) that the two bodies are far 
enough, so that they do not tidally interact; and ii) that the radiative energy loss 
does not affect the spins. Under these assumptions, the rates of the quadrupole 
tensor J, which allow to approximate the gravitational radiative energy loss 
(Eq. H22|) '). may be calculated with the system being simply modelled as two 
Newtonian point masses. Then we calculate that 



E : 



G if 2 



OPMjlim 



15 



2u2--(u.n)' 

4 ^ ' 



(30) 



where 



X = Xi — X2, u = X, K = Gmim2, r 



x/r 



(31) 



(mi and m2 are the zero-order masses of the pulsar and its companion; xi and 
X2 are their positions). Therefore, the Peters-Mathews coefficients are: 



ki — 2, k2 = - , fcdipoio = (scalar theory). 



(32) 
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In GR in harmonic coordinates, the expression of the energy rate is |34|: 

E ~ -(<Z'oPM)iim = -4 ^ ^ u2 - 11 (u.n)'l . (33) 

\b r"^ I J 

(I.e., the Peters-Mathews coefficients are in GR: 

fci — 12, fc2 = 11, fcdipoio = (GR in harmonic coordinates).) (34) 

In order to relate the energy rate to the rate of the orbital period P of the binary, 
one uses Kepler's third law, 

{2TT/Pfa^ = G(?7ii + 7712), (35) 

and the expression of the Newtonian energy of the binary, 

2a ^ ' 

(where a is the semi-major axis of the relative orbit), and gets 

------ (37) 

P 2E ^ ' 

For binary pulsars, usually the "companion" star is not seen and the observa- 
tional input is the list of the arrival times of the successive pulses emitted by 
the pulsar, in short "the pulse data" [H2!- Essentially, these data show a periodic 
variation due to the motion of the pulsar around the mass center of the binary 
system, but there is a slow drift in this variation, due to the decrease in the 
period P. If one wants to test a given theory of gravitation, then the "timing 
model" used to fit the pulse data should be consistently and entirely based on 
that theory. If we apply Eqs. l|Hn| and (p?7| to a given binary pulsar, in taking 
all orbital parameters as they are obtained from a timing model based on GR 
in harmonic coordinates, then we will obtain a period change of the same order 
of magnitude as that found with the timing model based on GR, but signifi- 
cantly different from it - because the energy rate (|30|l differs from that in GR 
(Eq. 133(1 ). while keeping the same order of magnitude for a given orbit and a 
given K . In view of the foregoing sentence about the consistency of the timing 
model and the theory, this does not prove that the scalar theory is not able to 
accurately fit the pulse data (which is indeed the true test). It is worth to recall 
that, in contrast with this, most other alternative theories give a very different 
energy rate from that predicted by GR, usually they even predict the wrong sign 



4 Conclusion 

A scalar bimetric theory with a preferred reference frame has been motivated 
and summarized. An asymptotic framework has been introduced for the "post- 
Minkowskian" (PM) approximation outside and inside the bodies as well, by con- 
sidering a family of initial-value problems for a perfect-fluid system; this seems 
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to be new and could be useful in other theories. In particular, it has been shown 
that this PM scheme is more general than the corresponding post-Newtonian 
(PN) scheme and that, therefore, a simple and reasonable justification can be 
given for the procedure usually adopted in studies on gravitational radiation. 
That is, it is allowable to take the energy rate as given by the PM approxima- 
tion and to say that it represents (approximately) the rate of the Newtonian 
energy. Using that PM scheme, the 0PM approximation of the energy rate for 
an isolated system has been calculated for the investigated scalar theory. This 
is indeed an energy loss, and it contains no dipole term (and no monopole term 
either). Thus, the radiative energy loss has the same structure in this theory as 
in GR, it has also the same order of magnitude. This is usually not the case 
for alternative theories, even for those that look much closer to GR. The author 
feels it plausible that this simple scalar theory might be able to fit accurately 
the pulse data of binary pulsars. 

Note 

The founding papers of special relativity, by Lorentz, Poincare, Einstein, and 
Minkowski, can be found online on the web. Links may be found at 
http : / / geo . hmg . inpg . f r/arminj on/INTRO . html. 
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